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Abstract 

Spatially extended stationary and traveling states in the strongly nonlinear 
regime of convection in layers of binary fluid mixtures heated from below are 
described by a few-mode-model. It is derived from the proper hydrodynamic 
balance equations including experimentally relevant boundary conditions with 
a non-standard Galerkin approximation that uses numerically obtained, sym- 
metry decomposed modes. Properties of the model are elucidated and com- 
pared with full numerical solutions of the field equations. 
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I. INTRODUCTION 



Convection in binary miscible fluids like, e. g., ethanol-water is a paradigmatic system 
for studying instabilities, bifurcations, complex spatiotemporal behaviour, and turbulence. 
Furthermore, it is sufficiently simple to allow for controlled experiments and the governing 
field equations are well known. So recently a lot of research activities [1-13] have been 
directed towards investigating the enormous variety of pattern forming behaviour in this 
system. The richness of spatio temporal phenomena in binary fluid mixtures stems from a 
feed-back loop between the fields of velocity, concentration, and temperature. Let us start 
with the velocity field: The convective flow is driven by the buoyancy force field which itself 
is determined by variations of the temperature and of the concentration field. The latter 
are on the one hand generated via the thermo diffusive Soret effect by temperature gradients 
and on the other hand they are reduced by concentration diffusion and by mixing through 
the convective flow. Since these changes influence the buoyancy which drives the flow the 
feed back loop is closed. 

In this article we concentrate on spatially extended convection states of straight parallel 
rolls that occur either as a horizontally traveling wave (TW) or in the form of stationary 
"overturning" convection (SOC) rolls. In particular, we derive a model which is able to 
describe also stable TWs in the strongly nonlinear regime and which goes beyond earlier 
approximations that are applicable only to the weakly nonlinear regime. We present and 
explore here the properties of a minimal Galerkin model that is based on a symmetry 
decomposition of numerically obtained fields. It reproduces the typical bifurcation scenario 
obtained from simulations and experiments. 

For fluid parameters typically realized in mixtures of water and about 10 wt.% ethanol, 
an oscillatory, subcritical onset of convection is observed. It is connected by an unstable 
TW branch with a saddle node bifurcation giving rise to stable, strongly nonlinear TW 
states. At a certain Rayleigh number, the phase velocity of these waves vanishes and the 
SOC branch of stable stationary states can be observed. Along the TW bifurcation branch 
which is shown in detail in Figs. 1 and 4 the concentration changes its structure from lateral 
homogeneity and vertically linear layering in the basic state over plateau-like distributions 
in fast TWs to boundary layer dominated, slowly traveling waves. The contrast between 
two adjacent TW rolls is strongly related to the phase velocity of the TW and it vanishes 
with this velocity. So, SOCs do not show such a concentration contrast. In SOCs adjacent 
rolls are mirror images of each other and they are separated from another and from the top 
and bottom plate only by thin boundary layers. The latter are a typical phenomenon for 
convection of weakly diffusing scalars. To capture this behaviour of the concentration field 
our model contains more degrees of freedom of the concentration field than of the velocity 
and of the temperature field. 

Our article is organized as follows: The second section presents the system, the fields 
needed for its description, their governing equations with the explanation of the relevant 
fluid and control parameters. The third section shows which symmetries in TW and SOC 
convection occur and how they change at the bifurcations. Furthermore, a discussion of 
the symmetry decomposed balance equations and of their numerically found solutions is 
given. Out of this data we extract in the fourth section a basis for a non-standard Galerkin 
approximation. Its results are presented in bifurcation diagrams of different order parameters 
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in comparison with the correct solutions of the system. We also discuss predictions of our 
model concerning the dependence of the bifurcation topology on the strength of the Soret 
coupling. 



II. SYSTEM 



The fields needed for the description of straight parallel rolls in a horizontal fluid layer 
exposed to a vertical gravitational acceleration g are temperature T, concentration C, and 
velocity u = (u,0,w). Here, u is the velocity field in the lateral x-direction and w is the 
velocity in the vertical z-direction. Fluid parameters are the Prandtl number a = u/k, the 
Lewis number L = D/k, and the separation ratio ip = — (3kx / '(cxT ) . v is the kinematic 
viscosity, k the thermal diffusivity, D the diffusion coefficient, a = —(l/p)(dp/dT) and 
(3 = — (1/p) (dp/dC) thermal and solutal expansion coefficients, p the fluid's density, kx the 
thermodiffusion ratio, and T (Co) the mean temperature (concentration) in the system. We 
scale lengths by the height d of the fluid layer, times by the vertical thermal diffusion time 
d 2 /n, temperatures by the imposed vertical temperature difference AT, and concentrations 
by AT a/ (3. Our control parameter is the Rayleigh number R = agd 3 AT/(uK,). Mostly we 
use the reduced Rayleigh number r = R/R® where R° c denotes the critical Rayleigh number 
for the onset of convection in a pure fluid (R® = 1707.762; R° c ~ 1747.28 in our model). 

We start with the hydrodynamic field equations in a nondimensionsal version using the 
scales defined above. They represent balances for momentum, heat, and concentration 



Assuming an incompressible fluid — this is reasonable as long as the convective velocities in 
the fluid are small in comparison with the velocity of the sound — the density is constant 
and the continuity equation as the balance equation for the total mass leads to 



In the equations (2.1), the temporal changes of a field are given by convective transport of 
that field on the left hand side of the equation. Apart from the advective rate of change 
the momentum (2.1a) changes by pressure gradients (— Vp) and dissipative effects (o"V 2 u). 
The gravitational acceleration acts as an external force (R [(T — T ) + (C — Co)] e z ) on the 
fluid. At that point, density fluctuations due to temperature and concentration fluctuations 
may not be ignored since they represent the driving of the system. 

The heat balance (2.1b) takes advective (u ■ VT) and diffusive (V 2 T) transport into ac- 
count. In the case of the concentration balance (2.1c), the cross coupling of the temperature 
to the concentration (— LipV 2 T) is a relevant effect. It is called thermodiffusion or Soret ef- 
fect [1,2]. For negative separation ratios ip, it acts in such a way that temperature gradients 
induce antiparallel concentration gradients; positive values of ip lead to parallel ones. If the 
Soret coupling ip vanishes, as it occurs in ethanol water mixtures of a concentration of and 
28 wt.%, the mixture behaves just like a pure fluid. 



d t u + u • Vu = - Vp + a {R [(T - T ) + (C - C )] e z + V 2 u 

d t T + u • VT = V 2 T 

d t C + u • VC = LV 2 C - L^V 2 T. 



} 



(2.1a) 

(2.1b) 
(2.1c) 



V-u = 
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From these remarks, the stability behaviour may be inferred: For negative separation 
ratios, the unstable layering of the density is reduced, for positive ones it is increased com- 
pared with the situation in a pure fluid. Thus, mixtures with ip < (ip > 0) have larger 
(smaller) critical Rayleigh numbers than the pure fluid. Below these threshold values, a 
motionless (u = 0) state with the laterally homogeneous conductive profiles 

Tcond(z) = T — z 
C cond (z) = C -ipz 

of temperature and concentration, respectively, is a stable solution of eqs. (2.1). Above 
threshold this solution becomes unstable to convection. It is convenient to use 

6 = T(x, z- 1) - T cond {z) (2.2a) 
c = C(x,z;t)-C cond (z) . (2.2b) 

as convective deviations of temperature and concentration from their conductive profiles for 
straight rolls in y-direction. 

In this paper we discuss convective roll solutions of (2.1) which have a lateral wave 
number of k = ir being close to the critical values around 3.116. These solutions are either 
stationary or they represent traveling waves. The latter are stationary in a frame comoving 
with the phase velocity v of the traveling wave relative to the laboratory frame. In the 
comoving frame {y = for states being stationary in the laboratory frame), we introduce 
the streamfunction $(x, z) which describes the incompressible twodimensional convective 
flow 

u(x, z) = (u(x, z), 0, w(x, z)) = (-d z , 0, d x )$(x, z) (2.3) 

in the x — z plane perpendicular to the roll axes. The velocity field u(x, z) (2.3) guarantees 
incompressibility, V • u = 0, by construction. 

The traveling and stationary states are solutions of the time independent field equations 

V 2 $) = a [Rd x (6 + c) + V 4 $] (2.4a) 
J($ } 0) = d x <Z> + V 2 6 (2.4b) 
J($, c) = ipd x § + LV 2 (c - tp9) (2.4c) 

in the comoving frame. They can directly be deduced from (2.1). Herein, we used the 
Jacobian 

J($, f) := (d x $)d z f - (d g $)d x f = (u ■ V)/ 

to describe the convective nonlinearity. 

The boundary conditions at z — ±1/2 in the comoving frame are for all x 

no slip u = —d z § = —v , 

w = d x $ = o , 

perfectly heat conducting 6 = , and (2.5) 
impermeable d z {c — ip6)= . 
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Numerical simulations [13] and a full scale Galerkin analysis with many trigonometric 
modes for the fields [15] have shown that the contribution from lateral Fourier modes higher 
than the first one in the velocity field is only in a range of 1% of the first one and that 
these higher velocity modes are unnecessary for reproducing the bifurcation topology and 
the peculiar structure of the concentration field, respectively. The analyses [13,15] have also 
proved the irrelevance of a lateral meanflow in the lateral velocity field. Therefore, we start 
our approximations with the ansatz 

z) = m ^ sin kxCi(z) + vz =: z) + vz (2.6) 

for the streamfunction. Here, C\{z) is the first even Chandrasekhar function [16], u> max 
the maximal flow amplitude, and v the phase velocity of the traveling state. These two 
velocities are the only degrees of freedom of the velocity field. By the above ansatz for 
the streamfunction we fix the origin x = 0, z = of the cooordinate system to a point of 
maximal upflow velocity. 



III. SYMMETRIES 

In this section we discuss symmetries of the fields in the coordinate system defined in the 
last section and we perform a symmetry decomposition of the fields and their balance equa- 
tions. After this, we analyse the relative weight of contributions from different symmetries 
in the solutions of the field equations. 



A. Decomposition 

Nonlinear states of stationary convection (SOCs) have definite parity [13] under the 
operation x — > — x: positive parity for the fields w, 8, and c, and a negative one for $ and 
u. Only the linear critical fields at the stationary bifurcation threshold have in addition 
definite vertical parity under the operation z — > —z: positive parity for the fields w, $, 9, 
and c, and a negative for u. 

Traveling wave states with non vanishing phase velocity have no definite lateral parity. 
Only that TW just on the merging point r* of the branches of SOCs and TWs has, of course, 
the same symmetry as the stationary states. At the point of the Hopf bifurcation out of the 
basic state, the linear critical fields have definite vertical parity, which is destroyed by the 
action of the convective nonlinearity in every nonlinear state. 

The temperature and concentration fields are decomposed in four parts each, correspond- 
ing to the four different combinations of vertical and lateral parity 

e = e ++ + e + - + e~ + + e— (3.ia) 

c = c ++ + c +- + c -+ + c — . (3.1b) 

Here, the first (second) superscript denotes the parity in lateral (vertical) direction. The 
symmetries themselves will be labelled by iS ±=l= in the same way. A symmetry decomposition 
of the velocity field or of the streamfunction is not necessary, since their spatial variation 
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is fixed by (2.4). The decomposition of the balance equation for the concentration (2.3c) 
yields the following four equations after inserting the streamfunction $ (2.4) into (2.3c) 



J($, c + -) - v d x c' + = LV 2 (c++ - V# ++ ) + *P d x $ (3.2a) 

J(Q, C++) - v d x c— = LV 2 (c+- - ^6+-) (3.2b) 

J($, c~) - v d x c ++ = LV 2 (c~+ - i)9- + ) (3.2c) 

J($,c- + ) -vd x c + ~ = LV 2 (c— -iIj6~~) . (3.2d) 

The decomposition of the heat balance equation (2.3b) can be performed in a completely 
analogous way. Using the ansatz (2.4) for the streamfunction the momentum balance can 
be split into two parts: one for the symmetry S ++ and one for S~ + 

-vd x V 2 $ = Rad x (V+ + (T+) (3.3a) 
= R d x (6 ++ + C++) + V 4 $ . (3.3b) 



All together, eight partial differential equations must be solved for the temperature and 
concentration fields. The two degrees of freedom of the velocity field, u> max and v, are 
determined by algebraic equations which result from projecting the equations (3.3) onto the 
fixed spatial modes of the streamfunction (2.4). 

B. Numerical solution 

The above equations have been solved by a multi-mode Galerkin method [15] that implies 
a full representation of the fields. In particular, it resolves the peculiar spatial structure 
in the concentration field [12,13] in a sufficient way. Fig. 1 shows the resulting bifurcation 
diagrams of TWs (full lines) and of SOCs (dashed lines). The TW solution branch bifurcates 
at the Hopf bifurcation threshold r osc with the Hopf frequency uu out of the quiescent heat 
conducting state and it merges at r* with the SOC branch. For these TWs we show in Fig. 2 
how the different symmetry decomposed contributions to the TW concentration field c(x, z) 
vary as one moves along the TW bifurcation branch w mgx = 0, i. e. from the bifurcation out 
of the conductive state, to the end of the TW states at r*. We measure the portion of a 
symmetry c ±± (x, z) in c(x,z) by (c++|c)/(c 2 ) using global averaging as the scalar product. 
Thus, c(x, z) contains at r osc (u> max = 0) only contributions of the symmetries S~ + (short 
dashed line) and S ++ (solid line). This reflects the fact that the critical eigenf unctions 
have positive vertical parity. However, already in weakly nonlinear states the portion of the 
symmetry S + ~ exceeds that of S ++ as a result of the strong increase of the zeroth lateral 
Fouriermode by the action of the convective nonlinearity. Near to the point where the TW 
phase velocity and the maximal convective velocity are equal and which has been identified 
[17] as the transition from weakly to strongly nonlinear TW convection the two symmetries 
S + ~ and S h have the same portion in c(x, z). Together, they make up 99.9% of it. Finally, 
at the SOC-TW-transition where the TW transforms into an SOC state there are no more 
contributions from S h and S . 

A correct stability analysis of the conductive state needs test functions of the symmetries 
S h and S ++ whereas the determination of the SOC-TW-transition which can be looked 
upon as a lateral parity breaking bifurcation of SOCs containing S ++ and S + ~ requires S~ + 
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and S . Therefore, in principle all symmetries are necessary to describe the two transitions 
in the system in an adequate manner. The yet unmentioned symmetry S (long dashed 
line in Fig. 2) contributes at most 0.06% of c(x, z) and is clearly the smallest portion in 
c(x, z) all over the bifurcation diagram. That is the reason why we think that modes of this 
symmetry — their influence will be discussed in Fig. 4 — are more or less irrelevant for the 
complex bifurcation scenario in the convection of binary fluid mixtures. 



The simplest approach is to consider only one mode per symmetry class and to discard 
the S part in the fields. This is the essence of older standard Galerkin approximations 
[14]. They yield linear relations between the square of the TW-frequency u, the square of 
the convection amplitude w max and the reduced Rayleigh number r — R/R® 



Here, ujh is the Hopf frequency at the oscillatory stability threshold r osc of the basic state. 
These linear relations reproduce the onset of oscillatory convection and its weakly nonlinear 
properties in an adequate manner. However, they are not able to describe a saddle node 
bifurcation and stable, strongly nonlinear TW convection. In fact, it was shown [18] that all 
Galerkin approximations using only one mode for each of the three symmetries S ++ , S~ + , 
and S + ~ yield linear relations like (4.1) — independent of the spatial shape of the modes. 

Because of the smallness of the portion of S in c(x, z) the simplest extension is to take 
two modes for each of the three relevant symmetries into account. This works reasonably 
well, as we will show below. Thus, such a model has to be looked upon as minimal model 
that is able to describe linear, weakly nonlinear, and strongly nonlinear states in binary fluid 
convection. 

In order to provide adequate modes we use numerical solutions of the system as a basis 
for our non-standard Galerkin approximation. As one is interested in the description of 
an upper stable TW branch and the saddle node bifurcation we symmetry decompose two 
states of the upper (numerically determined) TW branch and then orthonormalize them by 
a Gram-Schmidt-method using global averaging as scalar product. 

For fluid parameters that are typically realized in mixtures of ethanol and water (L = 
0.01, o~ — 10, ip — —0.25) we use a fast TW near the saddle (u = 4) and a slow one 
near the SOC-TW-transition (ou = 0.125). Both frequencies are small compared with 
the Hopf frequency of about 11.3. The resulting basis is displayed in Fig. 3. Since the 
structure of the temperature field changes only slightly along the whole bifurcation branch it 
is sufficient to consider only one temperature mode per symmetry. The three most important 
contributions to 9 are, like in the concentration field, the symmetries S ++ , S~ + , and S + ~ . 
Their contributions are very well described by the modes of the classical Lorenz model [19] 
and those of related approximations [14]. 



IV. NON STANDARD GALERKIN APPROXIMATION 



A. Mode selection 



J 1 = a(r - r osc ) + ^ 
wLx = P(r - r osc ) . 



(4.1a) 
(4.1b) 
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B. Nonlinear stable TW states 



Using the basis shown in Fig. 3 and the balance equations (2.3) one can derive a Galerkin 
model for convection in binary mixtures. This model has six or eight degrees of freedom 
for the description of the concentration field depending on whether the contribution of the 
symmetry S is contained or not. The temperature field is represented by three or four 
modes. As usual, projecting the field equations, e. g. (3.2), on the spatial modes of the 
respective fields yields a system of algebraic equations for the mode amplitudes. The order 
parameters of the velocity field, u> max and uj = vk, are computed from two algebraic equations 
for the velocity field that result after projection from (3.3). Here, we only want to discuss 
the bifurcation diagrams of the TWs, since those of the SOCs are well enough reproduced 
by earlier approximations except for quantitative values of M and their stability properties. 

Fig. 4 shows from top to bottom the TW bifurcation diagrams of Nusselt number 



These order parameters have been evaluated for extended models including S modes 
(dotted lines) and neglecting them (dashed lines) in comparison with the "true" solution of 
the system (solid lines). The first thing to notice is that for a very good reproduction of the 
stable TW branch (the upper one in N, the lower ones in M and ui) the symmetry S may 
not be neglected. Nevertheless, S is not necessary for a qualitative correct bifurcation 
scenario of a Hopf bifurcation, a saddle node bifurcation and a SOC-TW-transition on the 
stable, strongly nonlinear branch. This was not the case in earlier Galerkin models [14] 
connecting the basic state with the unstable SOC-branch by a TW-branch. 

The positions of the SOC-TW-transitions given by the model including the S sym- 
metry and neglecting it differ as can be seen in Fig. 4. This can be explained as follows. 
This transition may be interpreted as a stability threshold of SOCs whose lateral positive 
parity gets broken there. Since SOCs contain both S ++ and S + ~ they need to be tested 
for stability against disturbances of both S~ + and S . If one neglects S the stability 
analysis is incomplete leading to different values for the threshold. 

The correct bifurcation of the concentration variance M is due to an adequate description 
of the concentration field C(x, z) which is nearly in SOCs all over the convection cell except 
for small boundary layers between two adjacent rolls and at top and bottom plate so that the 
variance around its mean value is very small compared to the large concentration gradients 
in the basic state. 



The unstable TW branches resulting from the models differ from the exact ones [17], 
however without destroying the subcritical Hopf bifurcation topology — compare the dashed 



N = ((uT-VT)-e,) 



x\ z =-l/2 ~ 



1 - (d z 8) 



x\z=-l/2 ' 



TW frequency uj, and the reduced concentration variance 




C. Unstable TW solutions 
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and dotted lower unstable branches of N up to the saddle in Fig. 4a with the full line (for uo 
and M the upper branch corresponds to the unstable solution and is not shown in Fig. 4b 
and 4c). The reason for these differences is that the basic state is not tested for stability 
by the correct critical eigenfunctions. Since the oscillatory stability problem in binary fluid 
mixtures cannot be formulated as a variational priniciple as it is the case for the stationary 
one in a pure fluid [16] there is no need that an approximate stability threshold is higher 
than the exact one. 

As mentioned already the position of the Hopf bifurcation does not depend on the pres- 
ence of the symmetry S since the basic state needs only to be tested for stability against 
disturbances of the symmetries S ++ and S~ + . However, the initial slope of, e. g., the Nusselt 
number N being a nonlinear property involves modes of the symmetry S that are driven 
by the critical modes with the symmetries S ++ and S~ + via the convective nonlinearity. All 
this can be seen in Fig. 4a where the dotted and dashed lines have a common bifurcation 
threshold out of the basic state (N — 1) but spread in the nonlinear regime (N > 1). A 
qualitatively better reproduction of the unstable TW branch is obtained by changing the 
basis of the Galerkin model, namely by using states of the correct unstable TW branch. 
Then, the approximation quality of the stable states which are the main topic of this work 
deteriorates. 



D. Dependence on the separation ratio V 

Numerical simulations [13] have shown that the concentration distribution depends only 
slightly on the separation ratio ip for the small Lewis numbers L = 0(0.01) and large 
Prandtl numbers a = O(10) that are typically realized in mixtures of ethanol and water. 
This insensitiveness of the concentration field towards changes in ip suggests that our above 
presented models also describe convection in binary mixtures with values of tp different to 
those where the basis is extracted. 

The changes in the bifurcation topology caused by varying the strength ip of the Soret 
coupling may be described by tracing out the -^-dependences of the SOC and TW saddle 
node bifurcations r| oc and r^ w , and of the SOC-TW-transition r*. This is shown in Fig. 5 
for the model containing all symmetries (lines) in comparison with the full representation 
of the fields (symbols). Both saddles of the SOC and TW bifurcation diagrams procured 
by our non-standard model are in good quantitative agreement with the results of the full 
representation. The model only deviates for weak Soret couplings ip > —0.04. This can 
be explained by the smoother variations of the concentration field for those couplings (see 
e. g. [13]). The same holds for the SOC-TW-transition r* which in addition differs for 
strong Soret couplings ip < —0.3. This is due to the fact that the SOC-TW-transition can 
be interpreted as a boundary layer instability [20] which therefore is very sensitive to the 
thickness of the boundary layers in the concentration field. However, the thickness of the 
boundary layer is fixed in a model with fixed modes. Nevertheless, the model is able to 
describe the -^-dependence of both TW- and SOC-bifurcation scenario quite well. 
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V. CONCLUSION 

In this article we have presented a model describing traveling wave convection in binary 
fluid mixtures. Their complex bifurcation scenario is reproduced by a Galerkin method 
using two numerically obtained exact solutions of the system as a basis. These solutions for 
fixed fluid parameters are symmetry decomposed in order to meet the symmetry properties 
of the bifurcations and of the states in the system. We have shown that there exists a 
minimal number of degrees of freedom for the description of the concentration field. Our 
model uses this minimal number and reproduces numerically and experimentally observed 
states in binary fluid mixtures very well. In addition, it describes the dependence of the 
bifurcation topology on the strength of the Soret coupling quite well. Furthermore, it shows 
that in the temperature and velocity field the modes of the classical Lorenz model [19] are 
sufficient. Finally, the model gives useful hints for a derivation of an approximation using 
analytical test functions instead of the here applied numerical ones. 
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FIG. 1. Bifurcation diagrams of reduced convection flow intensity (a) and frequency (b) in TW 
(solid lines) and SOC (dashed lines) solutions in a binary mixture with tp = —0.25, L = 0.01 and 
a = 10. 
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Figure 2 

Strongly nonlinear convection in binary fluids 
St. Hollinger and M. Liicke, Z. Phys. B. 



FIG. 2. The deviations c(x, z) of the concentration field from of its conductive profile are decom- 
posed into four symmetry classes having definite parity under the operations x — > — x and z — > — z. 
The relative weight of the respective components c ±± in the complete field c(x, z) are plotted as 
a function of the maximal flow amplitude w nmx between the Hopf bifurcation (w ma _ x = 0) and the 
SOC-TW-transition (w max = w max (r*)). Parameters are L = 0.01, a = 10, and ip = —0.25. 
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Figure 3 

Strongly nonlinear convection in binary fluids 
St. Hollinger and M. Lucke, Z. Phys. B. 

FIG. 3. Orthonormalized basis for the description of the concentration field in the x-z-plane. 
The wave number is k = ir so that the lateral peridicity length is 2. The four different symmetries 
are denoted by iS ±=t at the side of the four blocks containing two basis elements each. The upper 
ones of each block are real states of the system (uj = 0.125) whereas the lower ones are linear 
combinations of two states resulting from the orthonormalization procedure. 
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FIG. 4. Bifurcation diagrams of Nusselt number N (a), TW frequency to (b) and concentration 
variance M (c) in the framework of Galerkin models with numerical basis including the symmetry 
S (dotted line) or discarding it (dashed line) in comparison with results of a many mode analysis 
(solid line). The filled squares show the two states which where used to construct a basis for the 
models. Parameters are L = 0.01, a = W, ip = —0.25, and k = it. 
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Figure 5 
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FIG. 5. Phase diagrams of the here presented non-standard Galerkin model (lines) in com- 
parison with the results of a full representation (symbols). The basis of the model was fixed at 
ip = —0.25 as indicated by the grey bar. The ^"dependence of the SOC-TW-transition r* (solid 
line, filled squares), the TW saddle node bifurcation r^ w (dotted line, filled circles), and the sad- 
dle node bifurcation of the SOC branch rg OC (dashed line, open squares) is shown in a double 
logarithmical plot. Parameters are L = 0.01 and a = 10. 
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